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Abstract. We present recent results on a systematic method to calculate transport coefficients for a meson 
gas (in particular, we analyze a pion gas) at low temperatures in the context of Chiral Perturbation 
Theory (ChPT). Our method is based on the study of Feynman diagrams taking into account collisions 
in the plasma by means of the non-zero particle width. This implies a modification of the standard ChPT 
power counting scheme. We discuss the importance of unitarity, which allows for an accurate high energy 
description of scattering amplitudes, generating dynamically the p(770) and /o(600) mesons. Our results 
are compatible with analyses of kinetic theory, both in the non-relativistic very low-T regime and near the 
transition. We show the behavior with temperature of the electrical and thermal conductivities as well as 
of the shear and bulk viscosities. We obtain that bulk viscosity is negligible against shear viscosity, except 
near the chiral phase transition where the conformal anomaly might induce larger bulk effects. Different 
asymptotic limits for transport coefficients, large- A^c scaling and some applications to heavy-ion collisions 
are studied. 
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1 Introduction 

The analysis of transport properties within the Heavy- Ion 
Collision program has become a very interesting topic, 
with many phenomenological and theoretical implications. 
Transport coefficients provide the response of the system 
to thermodynamic forces that take it out of equilibrium. 
In the linear approximation, energy and momentum trans- 
port is encoded in the viscosity coefficients (shear and 
bulk) whereas charge and heat conduction produce elec- 
trical and thermal conductivities, respectively. 

A prominent example of physical applications to col- 
lisions of heavy ions is found in viscosities. Although the 
matter produced after thermalization behaves as a nearly 
perfect fluid |T], there are measurable deviations, which 
are seen mainly in elliptic flow and can be reasonably ex- 
plained with a small shear viscosity over entropy density 
ratio [2] . In these analyses bulk viscosity is customarily ne- 
glected, based on several theoretical studies. However, it 
has been recently noted [3l|4] that the bulk viscosity might 
be larger than expected near the QCD phase transition, by 
the effect of the conformal anomaly. On the other hand, 
shear viscosity over entropy density is believed to have 
a minimum in that region. In that case, i.e., if the two 
viscosity coefficients are comparable at the temperatures 
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of interest, there are several physical consequences such as 
radial flow suppression, modifications of the hadronization 
mechanism or clustering at freeze-out [S]. 

Lattice analyses of transport coefficients are cumber- 
some, since they involve the zero-momentum and energy 
limit of spectral functions [6 , 7 , 8 , 9J , and these calculations 
are still not conclusive. It is therefore very interesting and 
useful to consider regimes accessible to theoretical analy- 
sis in order to provide complementary information about 
transport coefficients. The theoretical approach to trans- 
port coefficients has been traditionally carried out within 
two frameworks: kinetic theory and the diagrammatic ap- 
proach (Linear Response Theory). The kinetic theory ap- 
proach involves linearized Boltzmann-like equations and 
has been successfully applied in high temperature QCD 
fTUl and in the meson sector [TTJ[T2], while the diagram- 
matic method has been developed for high-T scalar and 
gauge theories [T3|fT4]. In both formalisms, it is crucial 
to include accurately the coUisional width, identifying the 
dominant scattering processes in the plasma. In the dia- 
grammatic framework, we have recently studied transport 
coefficients within Chiral Perturbation Theory [T5j 151IT71 
I18j . pertinent for describing the meson sector at low en- 
ergies and temperatures below the chiral phase transition 
[12] . Our analysis shows that in order to include properly 
the effects of the thermal width, the standard rules of Chi- 
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ral Perturbation Theory have to be modified for this type 
of calculation. 

In this work, we will present a detailed update of our 
formalism and main results, paying a more detailed at- 
tention to several aspects of formal and phenomenologi- 
cal interest. In particular, we provide a thorough deriva- 
tion of the relevant formulae, emphasizing the link be- 
tween Kubo's formulae and the definition of transport co- 
efficients through thermodynamic forces and fluxes. Our 
ChPT diagrammatic formalism is reviewed, describing the 
different contributions in the modified power counting and 
discussing specially the role played by the light resonances, 
which we generate dynamically by unitarizing with the In- 
verse Amplitude Method. Detailed results for all transport 
coefficients are given, underlying the connection with ex- 
isting phenomenological and theoretical calculations, in- 
cluding the predictions of non-relativistic kinetic theory, 
which we meet in the very low-T regime. We will also 
discuss the large- A'c behavior, which will provide an in- 
teresting check of some of the results obtained. 



2 Chiral Perturbation Theory 

In order to describe the dynamics of the light mesons (pi- 
ous, kaons and the eta) we will use Chiral Perturbation 
Theory, which is an effective field theory of QCD for the 
low-energy regime [50] • It is based on the spontaneous 
symmetry breaking of the chiral symmetry of the QCD 
lagrangian (for massless quarks): 

X = SU(3)l X SU(3)r = SU(3)v x SU(3)a SU(3)v , 

where the axial generators of the chiral algebra are broken 
leaving only the vector ones. As the result of this sym- 
metry breaking there must appear a number of massless 
Goldstone bosons equal to the number of broken gener- 
ators, and they are physically identified with the pions, 
kaons and eta. In order to construct an effective lagrangian 
it is necessary to obtain the transformation rules of the 
Goldstone bosons, 4>a, under the original chiral group. It 
can be shown [3D] that they transform non-linearly, so if 
we use the exponential parameterization for the Goldstone 
bosons we have: 

U{x) ^^ RU{x)L^ , with U{x) EE exp (^i*^^ , 
and (/)(x) = ^ (t)a{x)Xa , (1) 

a 

where R e SU{Nf), L e SU{Nf), Xa are proportional to 
the broken generators, and Fq coincides with the Gold- 
stone boson decay constant to lowest order in the chiral 
expansion (see below and [20j). This non-linear transfor- 
mation on U{x) implies the following transformations for 
the Goldstone bosons separately under the vector and ax- 
ial charges: 

[Qa, 06] = Uabc^c , [Qt 0b] = 9ab{<l>) , (2) 



where gabi<P) is some non- linear function. So we see that 
under the unbroken group the Goldstone bosons trans- 
form linearly but they do it non-linearly under the axial 
charges corresponding to the broken generators. Once one 
knows the transformation rules for the Goldstone bosons 
it is possible to construct a lagrangian which describes 
their dynamics as the most general expansion in terms of 
derivatives of the U {x) field that respects all the symme- 
tries of QCD: 

^chPT ~ C2 + + Ce + . . . , (3) 

where the subindex indicates the number of derivatives of 
the field U{x). In practice, we will deal only with £2 and 
£4, given explicitly by the expressions [TS] (for the Nf — 2 
case): 

£2 - ^ Tr{a^C/ta^t/ + Af2([/t + u)} , (4) 

and 

+ h2Tr{df,U^dM}T^r{d''U^d''U} 

- luM^ Tr{a^C/^9'^X/} Tr{U + C/^} 
8 

+ ^(^3 + k)M^{TY{U^ + U}f + h^M^ , (5) 

where Mg = 2Bom {ra = ~ is the quark mass) co- 
incides with the mass of the pion squared to lowest order. 
The coupling constants Fq, Bq, li and hi are called the 
low- energy constants, and are energy- and temperature- 
independent by construction. 

In order to deal with this infinite lagrangian we need 
a way of estimating the contribution from each Feynman 
diagram of interest, because we do not have an explicit 
coupling constant. Given a particular scattering ampli- 
tude M.{'mq,pi), where niq is the mass of the quarks (we 
will consider m = mu — rrid) and pi the meson external 
momenta, the dimension D of the diagram is defined by 
rescaling these parameters in the following way: 

M{tp,, t^ruq) = t^M{p.„ niq) . (6) 

Then, the dimension of a particular diagram can be easily 
computed (Weinberg's Theorem): 

D = 2 + ^Nn{n-2) + 2L , (7) 

n 

where L is the number of loops in the diagram, and Nn the 
number of vertices coming from the lagrangian £„. This 
dimension so defined actually tells us that the contribution 
from a given diagram is 0{{p/ A)^), where p represents 
an energy, momentum, meson mass or temperature. The 
scale A will be of order Ae ^ ^t^Ft^ ~ 1.2 GeV for energies. 
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momenta or meson masse^ and of order At ~ 300 MeV 
for temperature^ Therefore, the chiral expansion will be 
more reliable as we go down in energies and temperatures. 



3 Transport coefficients 

For a system out of equilibrium there exist thermody- 
namic forces (gradients of the temperature, the hydrody- 
namical velocity or the particle density) and fluxes, where 
the latter try to smooth-out the uniformities produced 
by the former in order to restore the equilibrium state 
of the system. Transport coefficients are defined as the 
coefficients for a series expansion of the fluxes in terms 
of the thermodynamic forces. Here we will deal with the 
transport coefficients corresponding to a linear expansion, 
specifically the shear and bulk viscosities, as well as the 
thermal and DC conductivities. According to relativistic 
fluid mechanics 22J, the energy-momentum tensor (flux 
of four-momentum) of a fluid can be decomposed into a 
reversible and an irreversible part: 



where 



(/; + hA''^N,)U''] 



(8) 



(9) 
(10) 



with e the energy density, P the hydrostatic pressure, 
the hydrodynamic velocity defined as the time-like four- 
vector that verifies U''{x)Uf,{x) = 1, Af'^ix) = g^"" - 
U^^{x)U^{x) is a projector, is the heat flow (difference 
between the energy flow and the flow of enthalpy carried 
by the particles), h is the enthalpy or heat function per 
particle, given hy h = {e + P) /n, N'^ is the conserved 
current (in case there is some in the system), and U'^'^ 
is called the viscous pressure tensor defined as the irre- 
versible part of the pressure tensor: 



The flows and U'^'^ can be written as: 



= {U^T"" - hN'')A^'„ , 



(11) 



(12) 
(13) 



^ It can be estimated as the momentum at which the contri- 
bution from a one-loop diagram with vertices of £2 that con- 
tributes to meson-meson scattering equals the tree level con- 
tribution from £,4 [5T]. It can also be estimated as the energy 
of the lowest resonance, Ae ~ 770 MeV, from where ChPT is 
expected to fail (see Section |6|. 

^ It can be estimated as the temperature corresponding to an 
average momentum equal to Ae, or as the critical temperature 
corresponding to the chiral phase transition. 



where we have split iT^"^ into a traceless part, 11'^'^, and a 
remainder (7J = — 7T'^^/3). If we now define the thermo- 
dynamical forces as: 



X 



vv^ — 



T 



hn 



^ (Zi^ A" + A" Af" )- ^A'^'A 

2^ cr T ' a r) ^ err 



(14) 
(15) 

r 

(16) 



with T being the temperature and V = A^'^d,y, then the 
shear viscosity, 77, bulk viscosity, C, and thermal conduc- 
tivity, K, are defined through the relations: 



n = cx , 



2r]X 



(17) 
(18) 
(19) 



At this point it is convenient to pass to the local rest 
frame (LRF) and recover the more familiar expressions 
for transport coefficients. In the LRF (we will denote by 
a tilde the quantities evaluated in this particular frame 
of reference), = (1,0), 77°" = (from the definition), 
and then we have: 



T, 



And for the case of the thermal conductivity: 

'aj; _ ckP^ 

T hn 

Now, using the relativistic Gibbs-Duhem relation 



T 



do 



hN' 



-kT 



n 



I 

T 



-TV ( - 
T 



we can rewrite (21 1 as 

f * - hN' 



(20) 
(21) 

(22) 

(23) 



Therefore we see that for a system without any conserved 
current (besides the energy-momentum tensor) the ther- 
mal conductivity is zero |23j . It is also convenient to ex- 
plicitly write the expression of the thermal conductivity 
for both the Landau's and Eckart's choices of the hydro- 
dynamical velocity [221 : 



Then, 



^N-{x)N,{x) 



rpOi 



(E) 



(L) 



~iQ T'^ 



(24) 
(25) 

(26) 
(27) 
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In the non-relativistic limit, we can neglect the second 
term in the right hand side of (21 1 and we recover, for the 



Eckart's choice, the Fourier's law T^-^^ — —ndiT. It is im- 
portant to remark here that these Landau's and Eckart's 
conventions apply to macroscopic averages of the currents 
T'^^ and over a fluid element, and not to the mi- 
croscopic currents themselves. The microscopic quantities 
will be relevant in the next section, where we obtain the 
expressions for transport coefficients in Linear Response 
Theory. 

Finally, for the DC conductivity, an electric current is 
induced in the gas by an external electric field which is 
constant in space and time, J* = a^-El^^. In general, the 
DC conductivity will be a tensor, but we will consider here 
the isotropic case, so cr^ = agij. 



3.1 Kubo's formulae for transport coefficients 

Let be a system described by a hamiltonian i/o (indepen- 
dent of time) to which we add a perturbation V{t), such 
that V{t) = for t < 0. Linear Response Theory (LRT) 
consists in taking into account only the linear effects pro- 
duced by the perturbation on the magnitudes of the sys- 
tem. For this to be a good approximation is necessary 
that V{t) be small in the sense that the eigenvalues of 
H{t) = Hq + V{t) and the eigenvalues i?!'^'' of Hq verify 
\Ea - £'i°V-^a ^ < 1- Then, it can be shown [24] that if 
0(t) is a certain operator in the Schrodinger's picture, the 
variation of the mean value of the operator produced when 
the perturbation is introduced is given (to linear order in 
the perturbation) by: 

6{d{t)) ^ {mmmt)) - {^^'\t)\d{t)\¥^\t)) 

D 

dt' {w('^\oMt-t')[6u{t),v^it'm^'\o)) , 

(28) 

where |tf'*-°^(0)) represents the state of the system at i = 0, 



|.f(o)(t)) 



»*|>f(o)(0)), \'P{t)) = e-'^*|!Z'(o)(0)), and 



On{t) = e'^«*0(t)e 



-'Hot 



(Heisenberg's picture). The re- 



sult ( |28[ ) is also valid if instead of mean values we deal 
with thermal averages ( • )t and then we want to study 
small deviations from thermal equilibrium. Note that in 
this case, according to (28), we calculate these deviations 



by evaluating the expectation value of the commutator at 
equilibrium. 

Now, by applying some small perturbation to our sys- 
tem in order to take it slightly out of equilibrium and 
using LRT, we can obtain the expressions for transport 
coefficients in terms of correlators. We start considering 
the DC conductivity. We perturb the system by coupling 
it to an external classical electromagnetic field: 



Therefore, the induced current is 



S{J,ix)) ^-ij d*x' G«U^ - 2;'Xxt(^') , (30) 

where G^^{x) = (?(i)([J^(a;), J^(0)]) is the retarded prop- 
agator. Now, using the gauge = 0, we get i?*^^ = 
— 9o^oxf And therefore, in momentum space, 

UJ ■' 

(31) 

Since the DC conductivity corresponds to the action of a 
constant electric field, and assuming spatial isotropy, we 
finally have: 



lim lim '-^^(^""^-^P) 
3 uj^o+ |p|->o+ UJ 

' lim lim ^^("'1^1) 



6 uj^o+ |p|^o+ oj 



(32) 



where p^r = 2Imi(G^)'j is the spectral function of the 
current-current correlator. The order in which the limit is 
taken is important, since the opposite would correspond to 
a field constant in time and slightly not constant in space, 
what would produce a rearrangement of the static charges 
giving a vanishing electric current. This is called a Kubo 
formula for the DC conductivity, and we can express it in 
another useful form in terms of the Wightman function 
G> by using the KMS relation G<{p) = e~'^'P°G>{p) and 
p=G> -G< mHH]: 

(T = -^ lim lim ^ I d^xeP-^[J\x),jm]) (33) 

= lim lim f d*x e'P-^ (r(x)MO)) , (34) 
3 w^o+ |pHo+ J 

where /3 = l/T, and we implicitly assume thermal aver- 
ages. As we show below, in the perturbative evaluation of 
transport coefficients the spatial momentum can be taken 
equal to zero from the beginning. 

Turning to viscosities, as we have seen, they are related 
to gradients of the hydrodynamical velocity in the fluid. 
Since we will evaluate the correlators at thermal equilib- 
rium, we can choose a global reference frame that is at rest 
with the fluid. We will give here a simple derivation of the 
Kubo formulae for the viscosities and the thermal conduc- 
tivity (for a more rigorous discussion see [25] ) . By perform- 
ing a boost that depends on the point, we can simulate 
gradients of the velocity, so the fluid velocity around some 
point xq becomes W{x) ~ x^ djU''{xo)- Then, this boost 
implies the change in the energy density STi — —U p 
where p is the density of momentum in the fluid. There- 
fore, this corresponds to the perturbation in the hamilto- 
nian density: 



V{t) = / Ux)A^^,{x) 



(29) 



V{x) = -T'^'x'djU 



(35) 
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Under this perturbation, the variation in the expectation 
value of the energy-momentum tensor is 

(36) 

where we have integrated by parts and used d^T^'^ — 0. If 
we now particularize for the case diJJ^ — 0, and compare 
with the expression (20), we then obtain: 



By choosing A^.^^ = 0, integrating by parts, and using 



In momentum space, we can write it as 
with 



1 lini lin, ^^>M 
20 w^o+ |pHo+ uj 



(37) 



(38) 



p,(a;, \p\) = / d^x e'P--([^.,(x),7r^^(0)]) , (39) 



and TTij = Tij — gijT^f^/i. In order to obtain the bulk 
viscosity we consider instead diUj — {l/3)dijdkU'' , and 
we have: 

5n = 5{V)=i J d^x' t' B{t - t'){[V{x),V{x')])dkU'' , 

(40) 



with V = —T\ /3. And comparing with (20 1, we get 



C = i / d^xte{t){[v{x),vm) . 



In momentum space we can express it as 

C=l lim lim P^^^ 
2 w->o+ |pHo+ lo 



with 



Pd^,\p\)^ d''xe^^-^{[V{x),Vm) 



(41) 



(42) 



(43) 



We now derive the Kubo expression for the thermal 
conductivity. As we have seen, it is necessary to have some 
conserved current in the system (besides T^^^) for the ther- 
mal conductivity to be non-zero. According to Eq. (23 1, 
heat conduction can be produced by a gradient in the 
chemical potential. In order to create such a gradient, we 
couple an external gauge field A'^^^ to the conserved cur- 
rent N^^ (^cxt pl^ys the role of an effective chemical po- 
tential), so the perturbation in the hamiltonian density 
is: 



(44) 



dfj_N'^ = we obtain: 



SP^S(f'^~hN'){x) 



(45) 



= i / dV t'd{t-t'){[r{x),N\x'mAl,{x') 



(46) 



with r EE T"* - hN\ Thus, by comparing with eq. ( |23| ), 
taking djA'^^^ constant, and assuming spatial isotropy, we 
have that the Kubo formula for k is: 



linr lim ^^^1^ , (47) 



6 ui^o+ |pHo+ lo 



p.iu, \p\) = / d^x e'P-mr{x),%m) ■ (48) 



with 



We here have used the following property of Wightman 
functions involving a conserved current N^^ (so d^N^^ = 
0): for any operator 0{x), and frequency w 7^ 0, it is 
verified 



d^x e''"*(iV°(a;)6(0)) =0 



(49) 



From the expression for we explicitly see that if there 
is no conserved current in the system (besides T^'^), then 
the thermal conductivity is zero. 

In other more rigorous derivations of the Kubo formu- 
las [25,26,27 , where they do not use energy-momentum 
conservation in the correlators, the expression of the bulk 
viscosity involves the operator 



/3 - c^T"" - pN" 



(50) 



where Cs is the speed of sound in the plasma and fi the 
chemical potential. Because of the property (49 1, this op- 



erator will give the same result for the bulk viscosity if it is 
calculated exactly (non-pert urbatively) . However, in our 
calculations we will use ( 50 1 and T' = — hN^ (instead 
of only —hN^) to obtain the bulk viscosity and the thermal 
conductivity, because we are interested in a perturbative 
calculation using propagators with a non-zero width, so 
the conservation of the energy-momentum tensor in cor- 
relators will not be exact within the level of appro xima tion 
we will use. As we will see, the extra term in (50 1 will 
be relevant in our approach. 



4 Transport coefficients in high-temperature 
quantum field theory 

At this point it is convenient to review what happens when 
one calculates transport coefficients in high-temperature 
field theories pUlfTB] . In these theories it turns out that in 
order to obtain the leading-order result for transport co- 
efficients a resummation of diagrams is necessary. As we 
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have seen, a transport coefficient is given in LRT by tak- 
ing the limit when the external momentum goes to zero of 
the imaginary part (spectral density) of some correlator. 
This process of taking the limit of zero external momen- 
tum implies the appearance of the so called pinching poles, 
which are products of retarded and advanced propagators 
sharing the same four-momentum: 



Gr(p)Ga(p) 



[Sip'' - Ep) + 5{p' + E,)] , (51) 



where Fp is the particle width (inverse of the collision time 
in the plasma) and Ep the particle energy. A pinching 
pole would correspond to the contribution from two lines 
in a diagram which share the same four-momentum when 
the external frequency is zero. For a theory at high 
temperature, F ~ ©(A^T), so ladder diagrams as the one 
depicted in Figure [l] all count the same, 0{1/X^), in the 
coupling constant and have to be resummed. 




Fig. 1. A ladder diagram with n rungs. 

Another kind of diagrams, bubble diagrams (Fig.|2]), in 
principle would give the dominant contribution, increas- 
ing with the number of bubbles according to the count- 
ing scheme given above, so they would be naively of or- 
der e'((l/A2)"A"-i) = e'(l/A"+i). But after some anal- 
ysis [T3] it can be shown that they can all be resummed 
giving a subdominant contribution (except for the bulk 
viscosity) with respect to the one-bubble diagram of Fig- 
ure [9| i.e., they correspond to the graph of Figure |3] with 
A-^ = + 0{X) and A2 = A^°\ with yl(o) being the 
lowest-order vertex. 

Fig. 2. A bubble diagram with n bubbles. 



Therefore it is interesting from the theoretical point 
of view to analyze what happens in ChPT, where we do 
not have and explicit coupling constant, in order to see 
whether a resummation is needed. 




Fig. 3. Generic representation of a ladder and/or bubble dia- 
gram in terms of two effective vertices, Ai and A2. 



5 Particle width 

As we have mentioned in the previous section, it is cru- 
cial to have lines dressed in the generic diagram of Fig. [3] 
(double lines) in order to take into account the collisions 
between the particles of the fluid, as dictated also by ki- 
netic theory. If the particle width was zero it would mean 
that particles would propagate without interaction, imply- 
ing that the corresponding transport coefficient would be 
infinite. We can approximate the interaction between the 
particles in the bath by considering the following spectral 
density with a non-zero width Fp-. 



p{p\p) = 



^ 

(po ~ Epf 



p2 

P 



(pO 



p2 

P 



(52) 



This approximation by a Lorentzian will be valid for a 
small enough width. The width is gencrically calculated 
for two-body collisions by [25] : 

^(^1) = '^™^4^^^^'* y di/2 dz/3 di^4 

X (2^)4^(4) (fc^ + k2-h- fc4)|T(s,t)|2 , (53) 

with di/, = d^k^/[{27Tf2Ei smh{l3Ei/2)], T{s, t) is the two- 
body scattering amplitude, and s = (fci + fc2)^ and t = 
(fci — fca)^ are the Mandelstam variables. 

If the gas is dilute, i.e. PE ^ 1 (dilute gas approxima- 
tion, DGA), the previous expression reduces to: 



r(fci) 



d^fc. 



(27r)3 



l_ P--E2/T 



a'tot(s)Wrol(l - Vi ■ V2) 



(54) 



where atot is the total pion-pion scattering cross section, 
the velocity of each the two colliding pions, and v^ei their 
relative speed. Up to energies of 1 GeV it can be shown [29] 
that for TTTT scattering only the channels / J = 00, 11, 20 of 
isospin-angular momentum are relevant and then we can 
make the approximation: 



O-tot(s) 



327r 
"3s" 

327r 



^(2J-Hl)(2/-|-l)|i/j(s) 



I, J 



[Koo(s)|2 + 9Kii(s)|2 + 5|i2o(s)|'] , (55) 



where t/j(s) are the partial waves, so the total scattering 
amplitude for tttt scattering is decomposed in terms of the 
isospin-projected scattering amplitude, Tj, and the partial 
waves as: 



T{s,t)^ 1^2(^1 +l)\Fj{s,t)\\ 



(56) 



with r/(s,<) = 327r^(2J+ l)i7j(s)Pj(cos6') , (57) 



j=o 



D. Fernandez-Fraile, A. Gomez Nicola: Transport coefficients and resonances for a meson gas in ChPT 



7 



and Pj being Legendre polynomials. Furthermore, in the 
00 and 11 channels there appear the /o(600) and p(770) 
resonances respectively. In order to deal properly with 
these resonances within ChPT we will have to unitarize 
our scattering amphtudes (see next section). The lead- 
ing order contribution to the pion width is represented 
by the diagram in Figure [4] Cutting this diagram in or- 
der to extract its imaginary part (the width) leads to the 
formula (53 1 with T{s,t) being the pion-pion scattering 
amplitude. All the diagramatic calculation will be carried 
out in the Imaginary Time Formalism (ITF) |24j which 
has the advantage of dealing with the same fields, vertices 
and diagrams as the corresponding vacuum field theory 
but with properly modified Feynman rules. 



e 



Fig. 4. Leading-order contribution to the pion widtli. 



6 Resonances 

It is difficult for ChPT to deal with the resonances that 
appear in some of the scattering channels because the uni- 
tarity condition 

SS^ ^l^lmtij{s)^a{s)\tij{s)\^ , 

with ct(s) = v/1 - 4M2/s , (58) 

is not respected for high-enough energy ^/s. This is be- 
cause the partial waves calculated in ChPT are essen- 
tially polynomials inp (and logarithms): i/j(s) = 4j (■?) + 



tfjis) + Cl(s^), with 47(s) = Ois") and s = 0{p'). In 
order to extend the range of applicability of the next- 
to-leading order results for the partial waves calculated 
in ChPT we will unitarize them by means of the Inverse 
Amplitude Method (lAM). The idea behind this method 
is essentially to construct an expression for the scattering 
amplitude which respects unitarity exactly and when ex- 
panded perturbatively matches to a given order the stan- 
dard ChPT expansion. The construction of this ampli- 
tude can be justified more formally by using dispersion 
relations [30ll3Tj . According to the lAM, the unitarized 
partial waves to order 0{p'^) are given by: 



Using this unitarization method, the /o(600) and p(770) 
resonances that appear in the pion-pion scattering chan- 
nels IJ = 00, 11 respectively are correctly reproduced 



for some set of values of the low-energy constants li (an 
overline denotes the renormalized low-energy constant, see 
([s])). In addition to appearing as peaks in the scattering 
cross section, resonances can also be identified as poles 
in the scattering amplitude after continued to the second 
Riemann sheet (SRS). Ift^^^ denotes the analytical contin- 
uation of the scattering amplitude off the real axis, then 
the scattering amplitude on the SRS, t'-^^\ is defined by 
ImiW(s-iO+) = ImtW(s + iO+), for s > AM^. Therefore 
one has 



1 -i2cr(s)t(i)(s) 



(60) 



A resonance corresponds to a pole of t'^^^ in the lower half 
complex plane, being the position of the pole related to 
the mass and width of the resonance by Spoic — {Mft — 
i//j/2)2, assuming that the resonance is a narrow Breit- 
Wigner one, which in the case of the /o(600) is not a so 
good approximation. Since we work in the center of mass 
reference frame, the mass and width obtained correspond 
to a resonance at rest. In what follows, we will fix the 
low-energy constants to the values Ti — —0.3, T2 — 5.6, 
= 3.4, Z4 = 4.3, which imply a mass and width at T = 
for the /o (600) /cr resonance of M„ = 441 MeV and 
Fa = 464 MeV respectively, and for the /o(770) resonance 
Mp = 756 MeV and = 151 MeV. 

If the pion gas is dilute enough so only intermediate 
two-pion states are relevant in the thermal bath, then we 
can define thermal scattering amplitudes as those calcu- 
lated like in the T = case but using thermal propagators 
instead The thermal scattering amplitudes i(s, T) 

(the temperature dependence starts to show up at 0{s'^)) 
also must verify the unitarity condition ( |58[ ) but making 
the replacement a{s) ^ aris.T) = cr(s)[r+ 2nB(-\/s/2)] 
(thermal phase space)rl there, and also in the expressions 
for the unitarized amplitudes. By considering these scat- 
tering amplitudes, we can study the evolution with tem- 
perature of the poles corresponding to the /o(600) and 
/9(770) resonances. In Fig.[5]we plot the evolution of the a 
pole with temperature, and we see that it remains broad 
at temperatures near the phase transition, while its mass 
is driven to the threshold as a possible indication of chiral 
symmetry restoration. 

Analogously, in Fig. |6] we plot the evolution of the p 
pole with temperature. We see that the mass shift is small 
remaining far from the threshold near the temperature of 
the chiral phase transition, and the width increases in- 
stead. This behavior is compatible with the dilepton data 
from heavy ion collisions (see details in f33, 341). 

We can also study the qualitative evolution of the res- 
onance poles in presence of nuclear density, introduced ef- 
fectively by varying the pion decay constant -FV, because 



(59) ^ {qq){p) 

Fm {mm 



1 - 0.35- 



(61) 



^ nB{x) = {e^'"^ — 1) ^ is the Bose- Einstein distribution func- 
tion. 
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T = 200 MeV 




T = 



250 MeV 



M (MeV) 

Fig. 5. Evolution of the /o(600)/(t with temperature in the 
mass- width plane. The red dashed line corresponds to the 
two-pion threshold. Each dot corresponds to an increment of 
10 MeV in temperature. The a remains broad near the phase 
transition although it experiences a big decrease in its mass. 



I\ (MeV) 



T = 2.50 MeV 




T = 200 MeV 



r = o 







M (MeV) 



Fig. 6. Evolution of the p(770) with temperature in the mass- 
width plane. The p increases its width with temperature and 
remains far from the threshold at temperatures near the phase 
transition. 



where p is the nuclear density, a-^N — 45 MeV is the pion- 
nucleon sigma term, and po ~ 0.17 fm~'^ is the saturation 
density of nuclear matter. This approach only takes into 
account a limited class of contributions, but we reproduce 
several aspects of the expected chiral symmetry restora- 
tion behavior at finite density. We have recently [34] im- 
proved the implementation of nuclear density effects for 
the /o(600) resonance by considering a microscopic cal- 
culation of many-body pion dynamics and unitarizing by 
solving the Bethe-Salpeter equation, obtaining in this way 
results qualitatively compatible with this simpler method 
we analyze here. In Fig. [7] we now plot the behavior of the 
tr pole at T = for several nuclear densities (the corre- 
sponding values of F^^ are indicated besides each pole) . We 
see that density effects drive the a pole faster towards the 
real axis, becoming a zero-width state, but the required 
nuclear density is very high (for instance, -FV = 55 MeV 



is equivalent to p ~ 1.9po)- At high enough density, when 
the pole crosses the threshold, there appear two separated 
poles on the real axis of the SRS below the threshold (vir- 
tual states), and for higher densities one of the two poles 
becomes a bound state (pole on the real axis in the first 
Riemann sheet). The virtual state which remains near the 
threshold and eventually becomes a bound state would 
correspond to a tttt molecule |34j . while the other virtual 
state behaves like the chiral partner of the pion in the 
sense that it tends to become degenerate in mass with it. 



-r/2 (MeV) 



bound state 




150 200 250 300 350 400 450 

M (MeV) 

Fig. 7. Evolution of the /o(600)/cr pole at T = for several 
nuclear densities controlled effectively by decreasing the pion 
decay constant. The corresponding values for F-n are indicated 
besides each pole. When the pole crosses the threshold it splits 
into two virtual states, and for higher densities one becomes a 
bound state. 



Analogously for the p(770), we see in Fig. [s] that den- 
sity effects also drive it faster toward the real axis and 
to the threshold. In this case however, after crossing the 
threshold, the pole becomes a pair virtual state-bound 
state located at almost the same value of mass and width, 
indicating a clear qq nature for this resonance |34j . 

In the following sections we will study the influence of 
unitarized scattering amplitudes on transport coefficients, 
as well as the influence of the in-medium evolution of res- 
onances on them. 



7 General analysis of diagrams for transport 
coefficients in ChPT 

In the analysis of transport coefficients within ChPT, anal- 
ogously to what happens in high-temperature quantum 
field theory, we also find non-perturbative contributions, 
oc 1/F (and 7^ = 0{p^)), due to the presence of pinch- 
ing poles. This would indicate that the standard ChPT 
power counting, dictated by Weinberg's formula ([7|, has 
to be modified in some way because otherwise, naively, di- 
agrams with a larger number of pinching poles would be- 
come more important as the temperature is lowered. We 
will show that for low temperatures, ladder diagrams are 
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Fig. 8. Evolution of the p(770) pole at T = for several 
nuclear densities controlled effectively by decreasing the pion 
decay constant. In this case, after crossing the threshold, the 
pole becomes a pair virtual state-bound state at almost the 
same position. 



the most relevant, but still perturbatively small in com- 
parison to the leading order given by the simple diagram 
of Fig. [21 




Fig. 9. Leading-order contribution to transport coefficients 
for low temperatures in ChPT. Double lines represent pion 
propagators with a non-zero thermal width. 



Again, the same topology arguments used in high- 
temperature theories are a priori applicable for the ChPT 
case, so we expect that the dominant contribution to trans- 
port coefficients come from ladder and bubble diagrams. 
We start by analyzing the spectral density corresponding 
to ladder diagrams. The spectral density of a generic dia- 
gram of the form shown in Fig. [3] can be easily calculated 
in ITF [T3| obtaining: 



lim Ptl^ = C(3 



*0+ CO 



2tt 



xGR,(e,p)GA(^,p)^i(f,p)/l2(C,P) , 

(62) 

where C is some combinatoric factor which depends on the 
kind of external field we couple to the pion loop, and it 
can be shown [15| that when considering a non-zero width 
we can take the zero spatial-momentum limit from the be- 
ginning. In the case of the simple diagram without rungs 
of Fig. [9| and for some constant external insertions, at 
T ^ we obtain that the spectral density behaves like 
liTa^^Q+ p{u>)/uj \/ Mt^/T, indicating that there could 
be important non-perturbative contributions from higher- 
order diagrams (ladder diagrams with an arbitrary num- 
ber of rungs) in the low-temperature regime. In order to 



give a first and naive estimation of the contribution at low 
temperatures (T ^ M„) from every diagram we assign a 
factor Y , that we expect to be of order Yi = ^M.,^lT for 
very low temperatures, to each pair of lines sharing the 
same four-momentum, and a factor X that we expect to 
be of order Xi = [A/7r/(47ri^7r)]^ for very low temperatures, 
to any other "ordinary" loop (Xi is the typical contribu- 
tion from a chiral loop). Therefore, according to this new 
counting, the contribution from a ladder diagram with n 
rungs would be of order 0{X'^Y^^^\ so ladder diagrams 
could in principle become more important as we go down 
in temperatures (where ChPT is expected to work bet- 
ter). Evidently, the contribution from the simple diagram 
in Fig. [9] would be of order 0{Y) instead of the 0(X) es- 
timation given by Weinberg's power counting. In order to 
verify this naive counting we have explicitly performed |15j 
the resummation of all the ladder diagrams for T <C 
and have found that it corresponds to multiply the lowest 
order result from Fig. [9] by some constant factor. This is 
because the contribution X from ordinary loops at very 
low temperatures is much smaller than naively expected, 
&oX ^Xi= ^TjM^ for T <^M^. Therefore the actual 
contribution from ladder diagrams at very low tempera- 
tures is 0(^2 Yi), so they are perturbatively suppressed 
when the number of rungs increases. However, although 
ladder diagrams give a contribution much smaller than 
naively expected, their actual contribution, ©(p^"), is still 
much larger than the estimated by Weinberg's counting, 
i.e., 0(p^"). Since in our ChPT lagrangian we also have 
derivative vertices, we expect that as the temperature in- 
creases, derivative vertices begin to dominate and the loop 
factor X increases over Xi . Also, the pinching pole factor 
would eventually become of order 0(1) and at temper- 
atures near the phase transition, in principle, we would 
have to resum all the ladder diagrams. 

Regarding bubble diagrams, it can be shown [T5^ that 
they all can be resummed giving a vanishing contribution 
in the limit w 0"*", for the leading order in l/F at very 
low temperatures. Other diagrams, like ladder or bubble 
diagrams with vertices coming from the lagrangians 
with n > 2, or diagrams with loops made with more than 
four-pions vertices would also be suppressed by the same 
arguments. But as we have commented before, as temper- 
ature increases Y becomes of order 0(1) and the particle 
width is not small anymore, so a resummation of all the 
possible diagrams would be in principle necessary. In this 
paper however, we will show the results corresponding to 
extrapolations of the leading order at low temperatures, 
and we will see that these extrapolations, when unitarized, 
give the correct order of magnitude for some observables 
near the phase transition, which indicates that the high 
temperature improvement due to unitarization is a key 
feature in this approach. 

We also remark that in the limit T ^ Af^, exponentials 
like that in (54 1 select small three-momenta of 0{\fM^). 



Therefore, in this limit it is enough to consider only the 
0{p^) amplitudes in the cross section, which allows to per- 
form a systematic T/M^r expansion fT5] for Fp and trans- 
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port coefficients, whose leading order results we give be- 
low. 

In the opposite limit we have massless pions (chiral 
limit) which is expected to be reached asymptotically at 
high temperatures. For = and using only the 0{p'^) 
amplitudes, the thermal width in (54) reduces to Fp = 
5T^|p|/(127r^i^^) and closed analytic expressions for trans- 
port coefficients can be given. However, these expressions 
should not be trusted at temperatures close to the phase 
transition, where three-momenta in the integrals are now 
of 0(T) and therefore the low-p expansion is not justi- 
fied and including unitarity corrections in the amplitude 
is crucial, as our results below will show. 



8 Electrical conductivity 

As we have seen, the DC conductivity is given in LRT by: 

lini lini P^^^ (63) 
6 w^o+ |pHo+ u} 

with 



obtained in lattice calculations [9 . The dots in the plot 
correspond to unitarizing the scattering amplitudes at fi- 
nite temperature, as we explained in Section [6] We see 
that the thermal evolution of the resonances does not af- 
fect much the conductivity. A more significant effect on 
the conductivity is produced when finite nuclear density 
is considered effectively by reducing F^^, but only at low 
temperatures. 



0.05 
0.04 
0.03 



p.(c., \p\) = / d^x e'P-%[J\x), Um ■ (64) 



Using the external sources method |20j . we couple an ex- 
ternal electromagnetic field to our ChPT lagrangian and, 
counting the electric charge e = 0{p/A), we calculate the 
lowest-order contribution to the DC conductivity, a^°\ 
given by the diagram of Fig. [g] Then, using the formula 
(62) we obtain [TS]: 



.(0) - 



d'^p IpI' 



ST J (27r)3 



E^F 



■nB{Ep)[l + n^{Ep)] , (65) 



where e is the charge of the electron and Ep = \J + 
For very low temperatures, T <C M^r, this expression adopts 
the simple form: 



r(0) 



15- 



e'Ft 



T1/2M, 



(66) 



It is interesting to compare our result with the expected 
kinetic theory (KT) behavior. According to KT |39j, tr ~ 
e^richT/M^ (rich is the density of charged particles, t is 
the collision mean time, and e is the particle charge), and 
T '-^ \ I F, F ^ nucr^Tr {"v is the mean speed of the parti- 
cles). In the non-relativistic limit, n ~ (VM^T)'^e~*^'"/-^, 
V ~ ^JtJM^, and constant, therefore a ~ 1/\/T. 

Thus, our result in ChPT is consistent with KT for T ^ 

In Fig. [To] we plot the lowest order contribution to the 
DC conductivity as a function of temperature for differ- 
ent choices of the scattering amplitudes that enter into 
the pion width. We see that unitarization (resonances) 
makes the conductivity grow from certain temperature. 
An increasing behavior for the DC conductivity is also 
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Fig. 10. Lowest-order contribution to the DC conductivity as 
a function of the temperature. The dashed line corresponds 
to considering non-unitarized partial waves to 0{p'^) in the 
pion width, the red line corresponds to the unitarized case to 
0(p*), and the dots correspond to unitarize at 0(p*) at finite 
temperature. The stars correspond to unitarized partial waves 
to 0{p'^) with Ftt = 50 MeV simulating nuclear density. 



As a phenomenological application of this result we can 
relate the electrical conductivity to the soft-photon spec- 
trum emitted by the gas of pions produced after a Heavy- 
Ion Collision (HIC) [TH]. The rate of photons emerging 
from a thermal system is related to the EM current-current 
correlator by: 



d'^p 



1 



^nB{uj)p^^{uJ= \p\) 



(67) 



Now, using the Ward identity Pf^p^'^ 
this rate to the conductivity as: 



0, we can relate 



di?^ 



1 



(--oMo|) = :i^3ra(T) 



(68) 



So the DC conductivity is directly related to the soft- 
photon spectrum, i.e., photons emitted with almost zero 
momentum. In order to compare with experimental re- 
sults we need to integrate this rate through the space-time 
evolution of the fireball produced after a HIC. For that 
purpose, we consider a simple hydrodynamical model of 
cylindrical symmetry in order to describe the expansion of 
the gas (Bjorken's model). Then, the measured rate would 
be given by: 



dN^ 
d'^p 



{PT -> + ) 



3r( r)a(r(r)) 
47r3 



rdr 



(69) 
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where we consider lead-lead collisions at SPS energies in 
order to compare with the WA98 experiment |35j , so -^/s = 
158A GeV, the nuclei radius is i?^ ^ 1-3 fm A^'^ ~ 7.7 fm, 
the expansion rapidity is ^r^nuci = 2 acosh(ys/(2A GeV)) 
~ 10.1, the initial proper time is ~ 3 fm/c, the final 
time Tf ~ 13 fm/c, and we consider the cooling law of an 
ideal gas T(r) = T,{T,/Tf/^ with T, ~ 170 MeV. With 
these parameters we obtain the estimate dN^ /d^p{pT — 
0) ~ 5.6 X 10^, that we indicate in the plot of the Fig. 
[TT] We see that other theoretical calculations do not fit 
the two lowest-energy points, while our result is compat- 
ible with a linear extrapolation to the origin from these 
two points. For low energies, the hadronic component is 
expected to dominate the spectrum, and in that regime 
a finite width in the particle propagator may be relevant 
due to the Landau-Pomeranchuk-Migdal effect, which was 
not taken into account in those studies and would make 
the contribution to the spectrum finite at the origin 156]. 



with = — hN'^ |^ But now, under this assumption, 
thermal averages would only imply sums over states of 
well-defined number of particles, |A^). Then, since in the 
diagram of Fig. [9] the energy-momentum enters to lowest 
order (i.e. without interaction), based on the KT theory 
expressions for the ideal-gas case [52]: 



d-^p 
(27r)3 

d^p 
(27r)3 



Epv'nBiEp) 



(72) 
(73) 



with — p'^/Ep, we define the operator iV^ through its 
Feynman rule for the vertex in momentum space heuris- 
tically as A^* = T'^^/Ep (the hmit of external momentum 
equal to zero is understood). According to this, the lowest- 
order contribution is then given by: 



, our estimation 



Correlation method 
Subtraction method 
^ Upper limit, 

subtraction method 




Pt (GeV/c) 



AO) 



1 



87r2r2 



d|p| 



nB{Ep)[l + nB{Ep) 



(74) 



As we have seen in the derivation of the Kubo formula for 
the thermal conductivity, here h represents the exact heat 
function per particle. However, we will approximate in our 
results h = (e + P)/n = Ts/n (s is the entropy density 
and n the density of particles) by the corresponding ideal 
gas expression, which we expect to be reasonable for the 
temperatures considered here. For very low temperatures, 
T < M^, we have: 



63 



^1/2^4 



(75) 



Fig. 11. Photon spectrum obtained by the experiment WA98 
[35j . We see that our estimate at the origin is compatible with 
a linear extrapolation from the two lowest-energy points. 



9 Thermal conductivity 

Although in the pion gas the only strictly conserved quan- 
tity is the energy- momentum tensor, in the energy and 
temperature regime we are dealing with, it is a good ap- 
proximation to assume that 2 —^ 2 collisions are the only 
relevant scattering processes, which in practice means that 
the pion number is approximately conserved [57] , yielding 
a nonzero thermal conductivity even when fi — [38'. 
Therefore, in order to compare with KT we apply 



, = lini lim P^^^ 
6 i^^o+ |pHo+ lj 



(70) 



where 



p,(u;, IpI) = / d'x &P-^[rix),%m) , (71) 



By KT [39], k ~ T^^{e — h)lv (e is the mean energy per 
particle and I ^ l/^aj^nfT') is the particle mean free path). 
In the non-relativistic limit, e ~ M^, h ~ 5T/2 -I- M^r, 
and then k ~ T^/^, so it is compatible with our result 



for low temperatures. In Fig. 12 we compare our results 



for K with a KT analysis [11 . Again, unitarity changes 
the behavior of the transport coefficient with temperature. 
We now see that density effects modify significantly the 
thermal conductivity at high temperatures, as expected 
when introducing and additional conserved charge (the 
baryon number). 



10 Shear viscosity 



It is given in LRT by: 
1 



20 



lim lim 

'^0+ |pHo+ 



Pni^^, \P[. 
to 



(76) 



* Note that in a previous calculation [17] we did not take 
into account the term with A'^', and therefore the result that 
we show here is qualitatively different at low temperatures. 
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Fig. 12. Lowest-order contribution to the thermal conductiv- 
ity as a function of the temperature. We compare with the 
analysis of [11], which is based on kinetic theory. 



with 

Pniu;, \p\) = J d\ e'P-^[r^{x),n,,m) > (77) 

and TT^^ = T*^ — 9^^T\ /3- Then, the lowest-order contri- 
bution is: 



(0) 



1 



d\p\ 



\P\ 



nB{Ep)[l + UBiEp)] , (78) 



For very low temperatures, T <ti Af^r, we have: 



V 



~ 37: 



(79) 



In Fig. [13] we compare our results for viscosities with those 
obtained by Prakash et al. using a KT analysis [IT]. We 
also agree with a work by Dobado et al. [T^] for the pion 
gas in KT. We see that nuclear density effects would only 
imply a significant change in the shear viscosity at low 
temperatures. By non-relativistic KT we expect the be- 
havior r?, C ~ Mt^vuI, thus ri,( Vt and both viscosities 
should then be of the same order at very low temperatures. 

Unitarity makes the quotient rj/s (s is the entropy den- 
sity) for the pion gas respect the bound lj(4n) predicted 



by Kovtun et al. [40], as we can see in Fig. 14 Without uni- 
tarization the Uncertainty Principle would be also violated 
eventually, since rj/s er/n ^ Et > 1. Furthermore, near 
Tc our value for rj/s is not far from recent lattice and model 
estimates [H]. Although we do not represent it in the fig- 
ure, we do obtain a behavior for rj/s growing slowly with 
T for temperatures (unrealistic) > 550 MeV. A slowly 
increasing behavior is also obtained by calculations from 
the quark gluon plasma (QGP) phase [12] , although a re- 
cent work predicts a more pronounced increase near the 
phase transition, in the so-called semi-QGP phase 
As another check, we can compute the sound attenuation 
length, which is given by (neglecting the contribution from 
the bulk viscosity) ~ 4ri/{3sT), and is directly related 
to phenomenological effects such as elliptic flow or HBT 
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Fig. 13. Lowest-order contribution to the shear viscosity as 
a function of temperature for several approximations in the 
scattering amplitudes. We also compare with the KT results 
of [H] (dots). 
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Fig. 14. Shear viscosity over the entropy density as a function 
of temperature. The horizontal dashed line corresponds to the 
AdS/CFT bound l/(47r). 



radn. We get, at T = 180 MeV, the value /; ~ 0.55 fm, 
in agreement with the estimate by Teaney |44j . 

In the chiral limit with only OijP'), we get exactly 
r] — 187rC(3)i^.^/(25r) with the Riemann's zeta function 
C(3) ~ 1.2. This 1/T decreasing behaviour, obtained for 
instance in [15], would imply that the AdS/CFT bound 
for 77/s is violated at some point, departing also from the 
phenomenological estimates discussed above. This high- 
lights the importance of reproducing correctly the high 
energy features of particle scattering, as we do within our 
unitarized approach. 



11 Bulk viscosity 

As we have seen, the expression for the bulk viscosity in 
LRT is: 



C^i lim lim ^^(^V^l) 



(80) 



ith 



Pc(/> \P\) = J d'x e'P-m'P{x),'Pm) , (81) 
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where, if there is no conserved charge in the system, V = 
~T\ /3 — c^Tf)Q. Then, the lowest order contribution to 
the bulk viscosity, corresponding to the diagram of 
Figure [9] is: 



d\p\ 



\p\H\p\y3-clE^^Y 



X nBiEp)[l + UBiEp)] , 

For very low temperatures, T <^ Af^, 
contribution simplifies: 



c 



(0) 



13- 



(82) 

the leading-order 
(83) 



In the chiral limit, — 0, we have the simple relation 
between the shear and bulk viscosities, = 15(1/3 — 
c^)^r]^^\ in agreement (parametrically) with the result ob- 
tained in the high-temperature regime of QCD This 
result implies that the bulk viscosity is suppressed with 
respect to the shear viscosity at large temperatures, as a 
consequence of conformal invariance, since = 1/3 for a 
free massless gas and, in fact, the bulk viscosity vanishes 
exactly for a conformally invariant theory. However, it also 
suggests that conformal breaking, as in the case of QCD 
through explicit and anomalous terms (see below) may 
induce sizable values for C/77, which would have interest- 
ing phenomenological consequences, since bulk viscosity 
is generically assumed to be negligible. This observation, 
supported by recent QCD analyses (see below) has led us 
to analyze recently in |18j the correlation between con- 
formal invariance and the bulk viscosity in the pion gas 
within the unitarized ChPT context. We reproduce some 
of the main results of that paper here, with more detail 
and emphasizing some quantitative aspects of the analy- 
sis. 



In the lowest-order contribution ( 82 1 , Cg is the exact 
(non-perturbative) speed of sound of the pion gas. How- 
ever, we can only calculate it within some approximation. 
In order to estimate the speed of sound we will analyze 
the relation between the trace anomaly and the bulk vis- 
cosity [SJHj. The scale invariance of the QCD lagrangian 
is broken explicitly by the quark mass and by the running 
of the strong coupling constant at the quantum level [?7] : 



25 



F^^^Fr + {I + l{g)}mqq , (84) 



where = T^'^Xi, is the dilation current, f]{g) is the /?- 
function, 7(17) is the anomalous dimension of the quark 
mass, and we consider the case of two flavors with m = 
rriu = rud- At finite temperature, the average of the trace 
anomaly is given in terms of thermodynamical quantities, 
{6)t = {T^^)t = e — 3P, and it has already been calcu- 
lated on the lattice for the pure glue theory [JS] as well as 
for QCD with almost physical quark masses [IH]. In Refs. 
[3] and |4] it was found through a low-energy theorem of 



QCD, and assuming some reasonable ansatz for the spec- 
tral function, a relation between the trace anomaly and 
the bulk viscosity. For instance, in the chiral limit, this 
relation reads [3]: 



do; 



P0e(t^,O) 



T 



d_ 
df 



{eh 



d (e-3P)* 



16|e.| 



(85) 



dT 

where (•)* = (•)t — {■)o is what is measured on the lat- 
tice, and \ey\ is the energy density of the vacuum. The 
spectral density peg involves the correlator of the trace 
of the energy-momentum tensor, and it is related to the 
pressure-pressure correlator by |71IB]: 



,{uj, 0) = 9ppp{Lo, 0) + uj5{Lu)TdT{e - 6F) 



7.0 



(86) 





K [Karsch ef a/., 03] 




1 ^HRG 
i 


il 


Lattice 




j ~ 170 MeV 






/ HRG with only tt, cr, p i 











6.0 
5.0 
4.0 
3.0 
2.0 
1.0 
0.0 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Fig. 15. Interaction measure calculated on the lattice (dots) 
and in the HRG approximation (green line, taking into account 
1026 states in total, see text), from [50|. The big red dot corre- 
sponds to the result from the HRG approximation taking into 
account only pions, the /o(600) and the p(770) states. 

Now, since the interaction measure, defined as Z\ = 
(6)* /T'^, has a peak near the critical temperature (see Fig. 
[l5| , it is reasonable to assume the following ansatz for the 
spectral density of the pressure-pressure correlator [HJHJ 
i: 



Ppp{lu,0) 



c 



UJUI, 







(87) 



where ~ 1 GeV. So a maximum in the interaction 
measure would imply a maximum in the bulk viscosity. 
Motivated by this, in order to obtain a good estimation of 
the speed of sound, we first calculate the interaction mea- 
sure in ChPT. The interaction measure can be calculated 
from the pressure through 

d 



dT 
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A calculation in ChPT of the pion gas pressure up to 3- 
loop order by Gerber and Leutwyler T9] is available. The 
diagrams contributing to the partition function of the pion 
gas are represented in Fig. 16 up to order 0(T^) according 
to the chiral counting. 



H O " oo o 



0(r») 



XCCOQ 

OO o o 



(X) 



Fig. 16. Diagrams contributing to order 0(T'^) and 0(r**) 
to the partition function of the pion gas. The first diagram 
just corresponds to a constant term in the lagrangian, boxes 
labeled with the number n correspond to counter-terms from 
the lagrangian 



Using the result for the pressure from in Fig. 17 
we plot the interaction measure for several orders in the 
pressure. The first peak corresponds to the explicit break- 
ing of scale symmetry by the quark mass and it has also 
been obtained in other works |51j . Interestingly, at order 
0{T^) there appears another peak near the phase transi- 
tion, that would correspond to the gluon condensate con- 
tribution. We then see that, in the chiral limit, interaction 
kicks in at 0{T^) (the 0{T^) diagrams would only imply 
a renormalization of the mass in this case). For AI-^ = 0, 
the trace anomaly has a simple expression at this order 



270 



In^ 
T 



(89) 



with Ap ~ 400 MeV for our choice of k given in Section |6j 
The result is very dependent on the choice of k, which m 
our case encode important non-perturbative information, 
as discussed in [18^. However, we see that the result of this 
perturbative calculation for the pion gas still is a factor 
10 smaller than the lattice and fuU-HRG results, see Fig. 

M 

The Hadron Resonance Gas (HRG) approximation con- 
siders a free (non-interacting) gas which consists of all the 
baryonic and mesonic states up to 2 GeV, 1026 in total 
|50j . Therefore, in this approximation the interaction mea- 
sure is given by: 



A 



(g"3P)- 

rp4 



1026 



E 



(e, - 3P,)* 



1026 oo 

2^E'? 

1=1 k=l 



i=l 
k+1 



rp4 

3 



k 



-Ki{kl3mi) , (90) 



where gi denotes the degeneracy of the state, 77 = ±1 
depending on whether it is a boson or a fermion respec- 
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Fig. 17. Invariant measure calculated perturbatively in ChPT 
for the pion gas. Each order of approximation corresponds to 
a definite set of diagrams in the partition function, see Fig. [16] 
The big red dot corresponds to the HRG approximation taking 
into account only the pion, the /o(600) and p(770) states (see 
text). 



tively, and Ki is the modified Bessel function of the sec- 
ond kind. In Fig. 18 from [i^, the lattice results for the 



interaction measure are compared with the HRG results 
for temperatures below the phase transition. We see that 
the HRG approximation fits better the lattice results near 
the maximum of the peak, while as the temperature de- 
creases they start to separate from each other. This might 
be due to the values of the quark masses taken in [49] , 
since the HRG approximation and ChPT should coincide 
for very low temperatures. It is important to remark that 
although the HRG approximation gives a value for the 
interaction measure compatible with lattice results near 
the peak, it is a monotonously increasing curve, so it does 
not have the shape of a peak, which the ChPT calcula- 
tion does have because it includes the interaction between 
the Goldstone bosons (in the HRG the non-zero interac- 
tion measure comes from the explicit breaking due to the 
non-zero mass of the states). 

We are also interested in the infiuence of the in- medium 
/o(600) and p(770) resonances on the trace anomaly and 
eventually on the bulk viscosity. In order to do it, we cal- 
culate the interaction measure in the Virial Gas Approxi- 
mation (VGA) which allows us to introduce the unitarized 
scattering amplitudes in the dilute gas regime. According 
to the VGA, to the lowest order in the interaction, the 
pressure of the gas is given by |53j : 



where 



(91) 



B 



(n) _ gtVi 



n+1 



27r2n 



dp p e 



2„-n/3(Ei-m,) 



(92) 
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[Cheng et al., 08] 
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Fig. 18. Comparison between the lattice results (dots and 
continuous line) with almost physical quark masses and the 
HRG approximation (dashed line) for temperatures below the 
phase transition, from |49| . 



B" 



27r3 



OO 

J dE E^Ki{(3E) (93) 



Y,{2I+l){2J+l)6%s{E) , 



According to this, we would not expect a big change in 
the bulk viscosity due to the in-medium evolution of res- 
onances. In fact, this is consistent with the idea [fS] that 
chiral restoration is not the main source for the conformal 
anomaly peak described above (unlike the vanishing tem- 
perature of the chiral condensate in ChPT, the position of 
this peak does not change in the chiral limit). Neverthe- 
less, a recent work based on the Linear Sigma Model, 
obtains a peak in the bulk viscosity at the chiral phase 
transition due to a minimum in the sigma mass, since in 
this model C oc ra/m"^- 

From the pressure we can also calculate other thermo- 
dynamical quantities like the entropy density, s = dPj dT, 
the specific heat, c„ = de/dT = Tds /dT , and the speed 
of sound, Cg = dP/de = s/cy. In Fig. 20 we plot the spe- 
cific heat and speed of sound squared for the pion gas in 
ChPT for several approximations. We see that the max- 
imum at the interaction measure implies a minimum in 
the speed of sound, so looking at the expression ( 82 1 we 



then expect a maximum in the bulk viscosity at the corre- 
sponding temperature. In Fig. [21] we also show the speed 
of sound squared and the equation of state as a function 
of the energy density obtained in the lattice with almost 
physical quark masses [49]. We see that the minimum in 
the speed of sound squared for the pion gas is still a fac- 
tor 2.5 less deep than the value from the lattice, but we 
have to bear in mind that we are dealing with a to 
two-flavor approximation, whose critical behavior should 
be that of a 0(4)-crossover. 



with $,i = Qpip-i-^i) ^ and S^jjg are the phase-shifts. In Fig. 
19 we plot the interaction measure for several scattering 
amplitudes in the VGA. The interaction peak is obtained 
when considering 0{p'^) amplitudes, and its height is al- 
most equal to the one obtained with the perturbative cal- 
culation of Fig. [l7| We then observe that the in-medium 
evolution of the /o(600) and p(770) resonances (see Sec- 
tion [6|) does not change significantly the height of the in- 
teraction peak. 




thermal unitarization, 0{p^) 



uiiitarized, 
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Fig. 19. Invariant measure in the VGA approximation (see 
text) for several scattering amplitudes. The in-medium evolu- 
tion of the /o(600) and p(770) resonances do no change signif- 
icantly the height of the interaction peak. 
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Fig. 20. Heat capacity (left) and speed of sound squared 
(right) for the pion gas in ChPT from the pressure obtained 
in [19]. At 0{T^), a maximum (minimum) in the heat capac- 
ity (speed of sound squared) is obtained near the chiral phase 
transition temperature. 



Finally, in Figs. 22 and 23 we plot the results in 



of the lowest-order contribution to the bulk viscosity and 
the bulk viscosity over the entropy density respectively. 
We explicitly show the importance of introducing unita- 
rized scattering amplitudes (resonances) in order to re- 
produce the peak near the phase transition. Nuclear den- 
sity effects do not change significantly the height of the 
anomalous peak, as expected from the previous analysis 
of the conformal anomaly. Comparing with Fig. |14[ our 
result for the ratio C*-°V'S is still smaller than rj'-^^/s near 
the transition, although the correlation with the confor- 
mal anomaly is clear, and that allows to predict larger 
C/s values if heavier states are included [18]. Neverthe- 
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Fig. 21. Lattice result with almost physical quark masses for 
the speed of sound squared and the equation of state, from 



less, these are extrapolations based on a low-T analysis 
and must be taken with a grain of salt. 



12 Large- ATc behavior of transport coefficients 

One of the main advantages of our formalism is that we 
can readily obtain the parametric dependence with the 
number of colors N^. This analysis is interesting, given 
the theoretical relevance of the large Nc to describe qual- 
itatively the QCD low-energy sector [55 . In addition, it 
will confirm some of our previous qualitative arguments. 

The large- counting of the low-energy constants li 
can be extracted from that of the SU{3) ones Li 
in the Nf = 2 hmit [56 , while = 0{Nc). This gives 
for the TTTT scattering amplitudes |rp ~ ©(l/iV^), regard- 
less of wheth er th ey are unitarized or not, and therefore, 
according to (53), we get Fp ~ 0{1/N^). This result, to- 
gether with the Nc scaling of the thermodynamic quanti- 
ties P~e~Cj~s~ C(l), which we extract from [T9] . 
implies that all transport coefficients scale as 0{N^) for 

However, in the chiral limit, from the expression of the 
pressure jT9] 
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Fig. 22. Lowest-order contribution for the bulk viscosity of a 
pion gas considering different approximations in the scattering 
amplitudes. 
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(94) 



and taking into account that log yip ^ li + 4^2 ^ 0{Nc), 
we get in this limit 1/3 - 0{lnAp/F^) - ©(l/TV,). 
Now, since (/rj ~ 0[{c^ — l/3)'^] in the chiral limit and 
the previous counting of the width is valid also in this 
limit, this means that for AI„ = 0, the scaling of the bulk 
viscosity is C ~ 0{1), unlike the other coefficients which 
still scale as O^N^)- Summarizing: 



For ^ 



(95) 



For = 



T^^^/s^O{Nl), C^C/5^0(1) • (96) 
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Fig. 23. Lowest order contribution for the bulk viscosity over 
the entropy density of a pion gas. The curves correspond to 
the same cases as in Fig. |22] 




These scaling relations are consistent with the results 
we obtained in the previous sections. The bulk viscosity is 
suppressed with respect to the shear viscosity in the chi- 
ral limit, as a consequence of scale invariance, although 
this is only a parametric dependence and it does not take 
into account the anomalous breaking near the transition 
[18]. For ^ the explicit breaking of conformal in- 
variance makes the two coefficients comparable, as we get 
for very low temperatures, where the mass terms domi- 
nate. For higher T, the chiral limit result is again reached 
asymptotically. Note also that the N^, scaling for Mj^ ^ 
is compatible with our leading expressions ( 66 ) , ([75| , ([79 1 
and 
for C 
Inyl 



(831 



above 



We disagree with the chiral limit TVc-counting 
given in |57j , where we believe that the scaling of 
discussed above is not properly accounted for. The 
Nc behavior is also useful in order to understand 



the origin of the different conformal-breaking terms near 
the transition [ISj. Finally, comparing with results from 
high-T QCD is also revealing. From the parametric ex- 
pressions given in [TU1H5] with the scaling as = 0{l/Nc), 
one gets ry/s ~ This is qualitatively compatible with 
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the idea of 77/5 approaching a minimum when coming from 
the low-T phase to the critical region, as we also obtain in 
our approach. In the high-T regime, C/*? is also suppressed 
by an additional (c^ — 1/3)^ factor. 



13 Conclusions 

Unitarized Chiral Perturbation Theory provides a con- 
sistent framework for the study of transport properties 
of meson matter. We have shown that, after a suitable 
modification of the standard ChPT power counting and 
including unitarity corrections in the scattering ampli- 
tudes in order to improve their high energy behavior, one 
ends up with a reasonable description of transport co- 
efficients for temperatures below the transition. At very 
low temperatures, our approach meets the predictions of 
non-relativistic kinetic theory, while at higher T we get 
an adequate behavior of transport coefficients when com- 
pared with existing studies based on the kinetic approach. 
In addition, we provide phenomenological predictions for 
the zero-energy photon spectrum and the shear viscosity 
to entropy ratio which are in fair agreement with data. To 
obtain these results, we have just considered the dominant 
diagram, with unitarized scattering in the thermal width 
for the internal pion lines. 

The results obtained within our approach for the bulk 
viscosity show a clear correlation with the scale anomaly, 
as suggested by previous works. Our formalism has the 
advantage of providing a theoretical analysis of transport 
coefficients for a massive pion gas without relying on lat- 
tice results, and therefore it might be useful in order to 
clarify the relation between the zero-energy limit of spec- 
tral functions involving the energy-momentum tensor and 
the bulk viscosity. 

We have also studied the large- A^c limit of the trans- 
port coefficients obtained in our approach. The paramet- 
ric scaling with Nc is consistent with our previous analysis 
and provides a qualitative description for the behavior of 
shear and bulk viscosities when approaching the critical 
region. 

Concerning future lines of research, we plan to in- 
troduce the effect of the strangeness sector (kaons and 
eta) which is relevant near the transition, where those 
states are no longer Boltzmann suppressed. The effect of 
pion chemical potentials, which has been sketched in our 
derivation here but not included in the results, is also an 
interesting extension of our work and will be considered 
elsewhere [37]. 
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